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ABSTRACT 


NASA is conducting a Geodetic Earth Orbiting Satellite-A (GEOS-A) 
Observation Systems Intercomparison Investigation. As part of this inves- 
tigation, some tests consisted of side-by-side tracking of the GEOS-A 
spacecraft by the Goddard Range and Range Rate (GRARR) system and the 
Goddard Laser tracking system. Seventeen passes were observed from 
July to November 1966 by the Rosman, North Carolina GRARR station and 
10 of these were evaluated. In the investigation, the Laser system tracks 
of the spacecraft were used as a reference trajectory for the GRARR system . 
The Laser data was smoothed using the GEOS-A Data Adjustment Program 
(GDAP) giving a reference orbit at the selected time of epoch in the form 
of a cartesian position and velocity vector. The evaluation of data shows 
that Laser orbits can be used to detect systematic errors in both range and 
range rate to about 2 meters and 1 centimeter per second respectively. 

Using the measured GRARR data and Laser orbital elements , GDAP 
determined the average range zero-set bias error to be -5.3 meters with 
a standard deviation of +2.5 meters per pass for seven of the ten passes. 
The other three passes were outliers with biases up to 30 meters from the 
established average. The range-timing error was determined to be 
-2 . 1 + 1 . 2 milliseconds . The unsmoothed range data RMS noise component 
was 6.8 meters after systematic errors were removed. It is assumed that 
inaccuracies in the transponder delay curve resulted in GRARR range-bias 
and timing errors. 

The range-rate residuals from GDAP were corrected for refraction 
and a sequential least squares regression program used to estimate co- 
efficients in various range-rate error models. No significant range-rate 
zero-set bias was detected. A consistent range-rate timing error of 
-0.20 _+ .02 milliseconds was observed, but its cause was not detected. A 
frequency-scale-factor error of about 10 parts per million was also ob- 
served. The unsmoothed range-rate data RMS component was 1 centimeter 
per second after these systematic errors were removed. 
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INTRODUCTION 


As part of the Geodetic Earth Orbiting Satellite-A (GEOS-A) Observation Sys- 
tems Intercomparison Investigation, several groups of tests were conducted from 
July 1966 through November 1966 at Rosman, North Carolina. These tests con- 
sisted of side-by-side tracking of the GEOS-A spacecraft by the Goddard Range and 
Range Rate (GRARR) system and the Goddard Laser tracking system. The tests 
were conducted to aid in the evaluation of the GRARR system and to determine the 
effectiveness of the Laser as a calibration instrument for electronic tracking sys- 
tems. A list of the spacecraft passes and a brief description of the data taken are 
given in Table 1. 


Table 1 

Spacecraft Passes 








Table 1 

Spacecraft Passes (Cont) 


Run No. 

Date 

Remarks 



anomaly was incorporated in the above analysis. For this 
type of data, if the two subgroups are combined after an 
ambiguity correction and a time correction are applied, the 
data may be useful . 

6 

Sept. 10 

This pass was not useful because of an unusual drift in the 
GRARR data . 

7 

Oct. 5 

No GRARR data was taken. 

8 

Oct. 6 

The Laser data suffered from a dropped bit equivalent to 

9 

Oct. 7 

thirty meters when certain data word configurations were 

10 

l 

Oct. 8 

encountered. This fault was quickly isolated to a gate mal- 
function in the range counter. The dropped bits were 
isolated by first smoothing all data, recognizing and de- 
leting the erroneous data points, and resmoothing the re- 
maining data. 

11 

Nov. 15 

This pass was quite short and tracking did not begin until 
after the point of zero range rate; hence, results were in- 
conclusive. 

12 

Nov. 18 

This pass showed an unusually large GRARR range bias. 

13 

Nov. 19 

This proved a good pass, but relatively short. 

14 

Nov. 19 

There was no overlap between the systems. 

15 

Nov. 20 

This was a good pass. 

16 

Nov. 20 

This was a good pass but had a larger than normal GRARR 
range bias. 

17 

Nov. 21 

This was a good pass. 


The Laser system tracks of GEOS-A were used as a reference trajectory for the 
GRARR system. The Laser data was smoothed using the GEOS-A Data Adjustment 
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Program (GDAP) 1 giving a reference orbit at the selected time of epoch in the form 

of a cartesian position and velocity vector. Then using the measured GRARR data 

and the Laser orbital elements, GDAP determined the range zero set error and the 

timing error for each pass . The range rate residuals from GDAP were corrected 

2 3 

for refraction and a sequential least squares regression program ’ used to esti- 
mate coefficients in various range and range rate error models. The results were 
analyzed to determine the appropriate models . 

LASER SYSTEM 

The Laser tracking system placed at Rosman for these tests uses an intense, 
highly-collimated, short-duration beam of light for illuminating the spacecraft 
being tracked. The beam is reflected at the spacecraft by special reflecting sur- 
faces and the returning light is detected photoelectrically, and its time of flight 
measured to yield range data. The actual laser transmitter was mounted on a 
radar pedestal along with a Cassegrainian telescope used for receiving the reflected 
Laser beam. When the Laser system is tracking, the transmitter is flashed at one 
pulse per second. Each transmitted pulse starts a time interval measuring unit 
necessary for range measurement. During the pass, the mount equipped with digi- 
tal encoders, is directed toward the expected position of the spacecraft by a pro- 
grammer fed with punched paper tape. By using a telescope, the operator can see 
the spacecraft and make corrections to keep it within the illuminating beam. Along 
with a range measurement, both the azimuth and elevation of the spacecraft are re- 
corded from the position of the mount. These angles are used for parallax correc- 
tions on the GDAP but not in the actual orbit determination since they need be accu- 
rate enough only for acquisition of the spacecraft and the driving of the mount during 
tracking. 

LASER SYSTEM PREPROCESSING 

The Laser tracking data consists of one measurement per second of range in 
nanoseconds and elevation and azimuth in degrees punched on paper tape. Using 
an AD/ECS-37 computer, punched cards are produced from the tape. Using a pre- 
processing program, the following corrections are added to the whole-second time 
for each measurement to give the time at the spacecraft. 
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• WWV correction (3.6 milliseconds for Rosman). 

• Delay time between the on-second pulse and actual Laser firing. 

• One-half the round trip interval of the Laser beam. 

The following correction is made to each range measurement when it is con- 
verted from a time interval into meters: 

• Internal delay correction due to photomultiplier, cables and receiving 
telescope (90 nanoseconds round trip). 

The corrected range, azimuth, and elevation measurements are put into a for- 
mat acceptable to GDAP. 

GRARR SYSTEM 

The GRARR system is a high-precision, spacecraft-tracking system that de- 
termines range using the sidetone ranging technique, and range rate applying the 
principles of coherent doppler. Angular data is obtained from X-Y mounted an- 
tennas but is not used for orbit determination. Each GRARR station uses an S-band 

4 

system and a VHF system in conjunction with a multichannel transponder on the 
spacecraft being tracked. Only the S-band system was used for this evaluation. 

The a-priori estimates of random errors are presented in Tables 2 and 3. 

GRARR SYSTEM PREPROCESSING 

Data at one measurement per second consisting of the range in meters, range 
rate in meters per second, and X and Y angles in degrees were used for this evalu- 
ation. In the operational preprocessing, the times are corrected to the spacecraft 
and a constant transponder bias correction equivalent to 3677 nanoseconds is made 
for each range measurement before submission to the Data Center at Goddard. For 
this evaluation the following corrections are made to each range measurement 
(R m ) in addition to the usual operational preprocessing described above: 

It should be noted that Table 2 is the same as Table 4 in the Goddard document 

g 

"Evaluation of Range Accuracy for the GRARR System at Rosman". However, the 
random range error due to transponder group delay variations has been omitted 
since this error is now modeled by the transponder delay curve used in the GRARR 
preprocessing. 
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Table 2 5 

GRARR Range Estimated Random Errors 


Name 

High Frequency 
(in meters) 

Low Frequency 
(in meters) 

Oscillator noise 

0.2 

0.0 

Thermal noise 

0.7 

0.0 

Quantization 

0.6 

0.0 

Digital timing 

4.5 

2.8 

Receiver delay variations 

0.0 

2.0 

Oscillator calibrations 

0.0 

0.5 

Transponder temperature 



and S/N variation delays 

0.7 

0.7 

Sub-total (RMS) 

4.7 

3.5 

Combined total (RMS) 

5. 

,8 


Table 3 

GRARR Range Rate Estimated Random Errors 


Name 

Random Errors 
(cm/sec) 

Thermal noise 

0.1 

Coherent oscillator instability 

0.1 

Oscillator noise 

1.0 

Quantizing noise 

0.5 

Coherent oscillator calibration 

0.0 

Digital timing errors 

0.4 

Total (RMS) 

1.2 
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• Additional uncorrected ground station system delays (Rosman, N. C. ) 

A R^ = 9.7 meters 

• Transponder delay (Channel A) 

AR 2 = (7.18 x 10 -8 )R 2 + (3.32 x 10~ 4 )R meters 

• Axis offset correction 

A Rg = 1.17 cos Y meters 

The final range (R T ) used is calculated by: 

R t = Rj^ + ARj - ARg + ARg meters 

No additional preprocessing corrections are made to the range rate data. 

There are two deterministic corrections which were not applied to the GRARR 
data at the time of processing. The first is the pre- and post-calibration adjust- 
ments which are now available. It is believed that these could affect some of the 
bias values by approximately one meter. The second is the WWV time synchroni- 
zation information. Because of the direct time link between the systems, this would 
not have affected the intercomparison results. 


GEOS DATA ADJUSTMENT PROGRAM (GDAP) 

After preprocessing has been completed, the GDAP determines a simultaneous 
solution to all the observation equations, which minimizes the weighted sum of 
squares of the residuals between these equations and the actual observations. Thus 
at any time T. we minimize: 

T ... -I 2 

F. = 7 w.f. (X., Y., Z., X., Y., Z., T., p, , . . . p ) 

J 4- L 1 * J J J 1 J J J 1 n J 


where f. = D., observed -D., corrected. The sum is taken over all data points. 

D. indicates the measurement, range, range rate or other data being processed. 
Xj, Yj, . . . etc. are the orbital parameters for position and velocity, T. = t^ - t Q 
where t Q is an epoch time, P^> • • • P n are error model parameters such as station 
location, zero set bias, etc. , and w. is a weight based on a-priori standard error 
inputs . 
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The implementation of the solution of the simultaneous equations depends upon 
the linearization of the equation set by expansion in a Taylor's series about a given 
approximation to the solutions. The resulting set of linear simultaneous equations 
is commonly known as the set of normal equations, and in matrix form is given by: 

NS = C 

where : 

N = B T WB 

and: T 

C = B WE 

The matrix B is called the matrizant and is the matrix of partial derivatives of 
the parameters at time tj with respect to the parameters at the time of epoch, t Q , 
and for the orbital parameters is represented by: 


B = 


ax, 


X 


az. 

X 


ax. 

z 

o 

dk 


The matrix W is the weighting matrix and is the inverse of the covariance ma- 
trix which is initially formed using a-priori input standard deviations of the param- 
eters. The vector E is given by: 

X., observed -X., corrected 
J J 

Z., observed - Z., corrected 
J J 

The normal equations are solved for the vector of parametric corrections, 8 
using matrix inversion. The solution, thus produced, is used to update or correct 
the current approximations to the error model coefficients . The solution is then 
iterated until the corrections are all equal to or less than one-half the a-posteriori 
standard error estimates, at which time the solution is complete. If the normal 


- 7 - 


equations are properly formulated, they yield as a by-product of the calculation, the 
inverse normal equation coefficient matrix which is given by: 


crX? <xX., Y. . . . <xX., Z. 
Ill 3 J 


V. = 
J 


<tY., X. <r Y? . . . o*Y., Z. 
J 3 3 3 3 


<xZ., X. crZ., Y. . . ,o-Z. 
3 3 3 3 3 


This matrix has the property of being the covariance matrix of the parameter 
correction values. This covariance matrix provides an indication of the remaining 
uncertainty in the determined values of the coefficients. 


LASER ERROR MODEL 

For the Laser passes, the GDAP state variables included the orbital elements 

• • • 

described by position and velocity X, Y, Z, X, Y, Z in cartesian coordinates at epoch 
time. The initial estimates of these parameters were obtained from Minitrack or- 
bital predictions and were unconstrained. Also, an azimuth bias and an elevation 
bias were used as state variables with an initial estimate of zero and an a-priori 
standard error of .5 milliradians for both. 


GRARR ERROR MODEL 

The orbital elements used for the GRARR passes were the converged values of 
the orbital elements from the same Laser pass. They were held constant by using 
a very small a-priori standard deviation. The radar error model included a zero- 
set bias parameter, initially estimated to be zero but allowed to vary freely by being 
given a very large a-priori standard deviation. There was also, in the error model 
a timing error term that was to detect any difference between the Laser and the 
GRARR range timing. This, too, was initially estimated to be zero since the dif- 
ference was expected to be no more than one hundred microseconds. Again this 
parameter was allowed to vary by giving the initial estimate a small weight. 


SUMMARY OF RESULTS 

Data were obtained from collocated Laser and GRARR systems at Rosman 
for seventeen passes. Tables 4 and 5 give a summary of the Laser and GRRAR 
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results, respectively. The average bias error for the Rosman GRARR relative to 
the Laser was found to be -5.3 + 2.5 meters. The random noise after removal of 
the bias was 6.8 meters. The average timing difference between the Laser and the 
GRARR range timing was -2.1 + 1.2 milliseconds. Later it is shown that the aver- 
age difference between the Laser and the GRARR range rate timing is only -0.20 + 
0.02 milliseconds. 

ANALYSIS OF ERRORS 

In order to investigate further the sources of error in the GRARR system, three 
passes were chosen for closer examination. The passes used were passes 10, 15 
and 17. Bias, RMS and range timing errors for these passes can be found in 
Table 5, while Table 6 gives average values for the three passes above. 


Table 4 

Summary of Laser Data 


Run No. 

Date 

Pass 

Duration 

(seconds) 

Range RMS 
(meters) 

Azimuth RMS 
(milliradians) 

Elevation RMS 
(milliradians) 

1 

July 27 

392 

27. 1* 

.34 

1.91 

2 

July 28 

252 

17.1* 

.17 

.11 

3 

July 29 

no overlap 

2.1 

.20 

.15 

4 

Aug. 10 

96 

4.3 

.14 

.09 

5 

Sept. 9 

345 

2.1 

.27 

.16 

6 

Sept. 10 

304 

1.5 

.20 

.07 

8 

Oct. 6 

205 

1.7 

1.07 

.14 

9 

Oct. 7 

377 

1.6 

1.11 

.24 


Oct. 8 

344 

1.2 

.26 

.13 

11 

Nov. 15 

174 

1.6 

.18 

.13 

12 

Nov. 18 

305 

1.0 

.31 

.07 

13 

Nov. 19 

153 

2.3 

.47 

.14 

14 

Nov. 19 

no overlap 

1.4 

.32 

.12 

15 

Nov. 20 

416 

1.6 

.69 

.30 

16 

Nov. 20 

336 

1.1 

.49 

.54 

17 

Nov. 21 

423 

1.6 

.51 

.24 

Average 


1.8 

.42 

.28 


*Not included in average. 
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Table 5 

Summary of GRARR Data 


Run No. 

Date 

Range RMS 
(meters) 

Range Bias 
(meters) 

Range Rate 
RMS 
(cm/sec) 

Range Time 
Difference 
(milliseconds) 

4 

Aug. 10 

12.9 

-5.6 + 1.0 

• 

No R data used 

-1.40 + .77 

8 

Oct. 6 

6.4 

8.7 + 1.3 

21.6 

-3.75 + .52 

9 

Oct. 7 

6.5 

10.6 + .8 

5.4 

-3.28 + .23 

10 

Oct. 8 

6.1 

-3.6 + .8 

3.3 

-4.04 + .28 

11 

Nov. 15 

5.6 

-4.1 + 2.2 

3.3 

-1.73 + 1.27 

12 

Nov. 18 

6.1 

-35.2 + 1.0 

8.6 

-.82 + .42 

13 

Nov. 19 

5.7 

-4.7 + 1.5 

9.9 

-.77 + .11 

15 

Nov. 20 

6.2 

-2.6 + .6 

2.0 

-1.47 + .17 

16 

Nov. 20 

6.2 

-10.1 + 1.1 

5.8 

-2.02 + .44 

17 

Nov. 21 

6.2 

-6.5 + .7 

2.4 

-1.41 + .17 

Average 

6.8 

-5.3 + 12.4 

6.9 

-2.07 + 1.19 


Range and range rate data is reported only for Laser passes considered accept- 
able. GRARR runs were made for all Laser passes where data was available; how- 
ever, as expected, results were poor where Laser or GRARR malfunctions were re- 
ported . 

It was felt that the standard deviation from GDAP associated with the bias fig- 
ures given in Table 5 do not represent the true standard deviation that can be ex- 
pected on any one pass using GRARR data with Laser defined orbits. Omitting the 
runs (8, 9, 12) with a bias which deviate greatly from the norm, a new average 
bias was calculated for the remaining seven passes. Coincidentally the new bias 
was calculated to be -5.3 meters. 

Using the formula: 
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where N = 7, B = -5. 3 meters and is the bias for the i pass. A more 

realistic standard deviation of + 2. 5 meters was calculated and will be reported as 
the per pass standard deviation of the range bias. 


Table 6 

Selected Data Set Averages 



Range RMS 
(meters) 

Range Bias 
(meters) 

Range Rate RMS 
(cm/ sec) 

Range Time Difference 
(milliseconds) 

Average 

6.2 

-4.2 + 2.0 

2.6 

-2.31 + 1.50 


When the adjustment program converged for a given pass, the range and range 
rate residuals obtained using the final parameter set were punched out on data 
cards. Since range rate refraction correction was not made in the GDAP at the 
time of these reductions, this set of range residuals was refraction corrected and 
a new set of range rate residuals obtained. The final refraction corrected range 
and range rate residuals was input to a sequential least squares regression analysis 
program which was used to determine the statistical feasibility of a proposed range 
or range rate error model. 

LASER ERRORS 

For the three Laser passes used for GRARR error modeling, the Laser range 
residuals had a mean value of -0 . 002 meters and a standard deviation of 1. 5 meters. 
A Chi-square test of normality was run on the range residuals and none of the three 
sets were found to be significantly non-normal, although a slight skewness was 
noticed in each of the three data sets. However, the residuals appeared to be ran- 
dom and if any systematic effects were present, they were quite small. 

7 8 

Tests have been made for serial correlation in the Laser data. ’ The results 
showed the serial correlation to be insignificant. Using the assumption of indepen- 

9 

dence, the eigenvalues of the covariance matrix of the velocity parameters X, Yand 
Z can be used to get an estimate of the accuracy with which range rate can be de- 
termined using orbits defined by the Laser. It was calculated that range rate would 
be no more accurate than one centimeter per second using the Laser as a standard. 
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TEMPERATURE AND S/N VARIATION OF GRARR TRANSPONDER DELAY 


Temperature and S/N variables in transponder delay were measured before its 
installation in the satellite. Table 7 summarizes these errors which may be re- 
moved by a look-up table as a function of temperature and up-link signal strength. 
These errors have not been removed from the data presented in this document. 


Table 7 

Transponder Variations with Temperature and Signal Level 


Run No. 

Date 

Temperature 
(degrees centigrade) 

Power Up-link 
1500/3000 km 

Range Variation 
(meters) 

4 

Aug. 10 

4.4 

-59/ -65 dbm 

+ .75 

8 

Oct. 6 

16.7 

-59/ -65 dbm 

+ .45 

9 

Oct. 7 

3.6 

-59/ -65 dbm 

+ .78 

10 

Oct. 8 

4.2 

-59/ -65 dbm 

+ .75 

11 

Nov. 15 

-4.6 

-59/ -65 dbm 

+ .88 

12 

Nov. 18 

5.1 

-59/ -65 dbm 

+ .73 

13 

Nov. 19 

-8.4 

-59/ -65 dbm 

+ 1 . 00 

15 

Nov. 20 

-7.9 

-59/ -65 dbm 

+ .98 

16 

Nov. 20 

4.9 

-59/ -65 dbm 

+ .74 

17 

Nov. 21 

-6.1 

-59/ -65 dbm 

+ .95 


Test specifications require that the maximum tolerance for transponder delay 
shall not exceed 32 nanoseconds or 4.8 meters over a signal level of -15 dbm to 
-85 dbm with the temperature between -10° C and + 45° C. 

GRARR ANTENNA POSITION BIAS 

The Rosman S-band GRARR parabolic antennas are mounted on an X-Y sys- 
tem. The X axis is oriented North-South and is fixed in position. The Y axis is 
oriented East-West, but is mounted 1.17 meters above the X axis. This causes 
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the position of the Y axis to vary sinusoidally in respect to the geodetic location of 
the site. The correction in meters is: 

A R = 1.17 (Cos Y), where Y is the Y angle of the pedestal. 

ROSMAN ZERO SET BIAS 

The average zero set bias for these tests after removal of known biases was 
-5.3 meters with a standard error of the bias estimate at 2.5 meters. Known 
biases, which were removed, include the 14-foot RE-142 cable used in precalibra- 
tion adjustments (at the collimation tower) which accounts for 6.3 meters. Also, 
the displacement of the horn aperture on the collimation tower reflector, consider- 
ing reflections within the cassegrainian feed system, is 1.0 meter greater than that 
previously used. The delay in the collimation tower transponder, when zero setting 
the system, has been found to be 3.0 meters greater than had been thought with re- 
spect to the test set used to calibrate both the collimation tower transponder and the 
satellite transponder. The 156 meters distance between the collimation tower and 
the tracking system antenna was measured to the stationary center of the X-axis of 
the GRARR antenna. The correction described in the preceding paragraph now 
accounts for the offset between the Y- and the X-axes , treating the data as if it 
were initially referenced to the center of the Y-axis and correcting it to the center 
of the X-axis. Therefore, the collimation tower data must be treated as if it were 
initially referenced to the center of the Y-axis . This requires an adjustment of 
0. 70 meters. These errors in the GRARR system at Rosman, North Carolina are 
listed in Table 8. 


Table 8 

GRARR System Ranging Errors for Rosman, North Carolina 


Name 

Range Error 
(meters) 

Zero set bias to X-axis (already corrected at Rosman 
prior to transmission) 

156.00 

Aperture of transponder horn 

-1.05 

Cabling delays from aperture to transponder 

-6.30 

Pole beacon transponder delay 

-3.00 

Boresight calibration conversion, X to Y axis 

+ 0.70 

Boresight parallax error due to dish separation 

-0.02 

Total uncorrected error 

-9.67 
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The net 9. 7 meters correction was made in the GRARR preprocessing pro- 
gram. 

TIMING ERRORS 

The average range timing difference between the GRARR system and the Laser 
was -2.1 + 1.2 milliseconds. However, this timing error appeared to drift since 
the errors for the October passes varied from -3.28 milliseconds to -4.04 milli- 
seconds with an average value of -3.7 + 0.4 milliseconds while the errors for passes 
in November varied from -.77 milliseconds to -2.02 milliseconds with an average 
value of -1.4 + 0.5 milliseconds. Also, using techniques of linear regression, a 
consistent range rate timing error of -0.2 milliseconds appeared as the coefficient 
of the range acceleration term in the range rate model. 

Part of the timing errors can be explained by two known lag delays in the cir- 
cuitry of the GRARR System. Both the range and the range rate are affected by a 
0.03 millisecond delay in the GRARR start pulse. Also, a narrowband, single- 
tuned LC circuit filter precedes the range discrimination circuit in the ground re- 
ceiver. The time delay through such a circuit may be expressed as: 



where Q is the quality factor of the tuned circuit (measured to be equal to 84) and 
f is the major side tone frequency (100 kc). 

The delay caused by this circuit is 0.13 milliseconds, giving a total explained 
delay of 0.16 milliseconds. This time delay has not been removed from the timing 
figures given in the rest of the paper. 

It should be noted that the difference between the GRARR range and the GRARR 
range rate timing errors seems to indicate that the timing problem was not caused 
by the Laser, but probably by the GRARR circuitry affecting only the range measure- 
ments . 

TRANSPONDER DELAY VS. DOPPLER DELAY 

The up-link receiver contains a filter which produces a phase delay according 
to the frequency of the input. All GRARR tracks of GEOS- A used in this experiment 
were taken on channel A at an input frequency of 2271.9328 megacycles. This is 
significant since each transponder channel has a characteristic delay curve. 
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In September of 1964, the Military Electronics Division of Motorola, Inc. tested 
the delay in the GEOS-A transponder and found it to be 3511 11 nanoseconds at zero 
doppler for channel A. The transponder was tested again at Cape Kennedy, Florida 
in November 1965 prior to the launch of GEOS-A and the delay was found to be 3677 
nanoseconds at zero doppler. When the GRARR data is received at GSFC, a cor- 
rection of 3677 nanoseconds (551.6 meters) is made to the range data taken on chan- 
nel A. The solid curve in Figure 1 shows a graph of the range error vs. range rate 
(doppler) based on the Motorola values. However, the entire curve has been verti- 
cally offset by a value of 166 nanoseconds (24.9 meters) to bring the delay at zero 
doppler to the 3677 nanosecond value used at Goddard. It was assumed that the 
shape of the transponder correction curve was not altered. A second degree poly- 
nomial was fit to the solid curve in Figure 1 to arrive at the transponder delay 
correction used in the GDAP preprocessor. The standard error of the fit for chan- 
nel A was + 0.2 meters over the operating region of GEOS-A which is between 
+ 50 kc doppler. 

Inasmuch as there seemed to be some uncertainty about the transponder delay 
correction, the range residuals were subjected to a regression analysis using a 
second degree polynomial of the form: 

(a) A R = K g R 2 + KgR 

Since this model is of the same form as the transponder correction, it was hoped 
that any error in the present correction could be resolved. The reduction in the 
rms value of the range residuals, however, was insignificant as shown in Figure 3. 

Also, it was felt that any transponder errors affecting the range measurements 
would affect the range rate residuals in a manner represented by a differential form 
of the range curve (a) given by: 

AR = 2K 4 RR + K 2 R 

• • 

where R is range acceleration. Here it was expected that the regression coefficients 

. 2 

and K 2 would be significantly similar to the coefficients of R andR, respectively, 
in the range transponder correction (ARg, page 6). Since this did not occur, a fur- 
ther search into the GRARR hardware characteristics was initiated. 
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Two oscillators are used in the GRARR ground station. One for the up-link 
carrier frequency used to extract doppler information, and the other for the rang- 
ing side tones which are phase modulated on the up-link carrier. Since the space- 
craft transponder delay has a different effect on signals traveling with group velocity, 

12 

such as the ranging side tones, and signals with phase velocity, such as the carrier , 
dissimilar delay characteristics of the range and the range rate are to be expected. 
The delay which affects the ranging side tones has been modeled; however, work is 
still in progress on the delay characteristics of the carrier through the transponder. 

The offset oscillator on the spacecraft has an accuracy specification (see refer- 
ence 10) of five parts in 10 or about 8. 5 kHz at 1705 MHz. Any drift in the oscilla- 
tor would alter the transponder curve given in Figure 1. This change in the trans- 
ponder delay curve could cause an inaccurate transponder delay correction in the 
spacecraft range calculation causing what would appear to be a timing error and a 
bias error in the orbital fit. The range bias errors caused by a frequency offset 
would be quite small and well within the standard error of the range bias noted in 
this paper. However, an inaccurate zero setting of the transponder delay curve in 
Figure 1 could introduce a significiant range bias into the orbital fit. For example, 
the range bias of -5.3 meters given in this paper could have been caused by a zero 
set error of approximately 35 nanoseconds in the transponder delay curve, i.e. , 
using, as a transponder delay correction, the solid curve in Figure 1, with a zero 
doppler correction of 3642 nanoseconds rather than the 3677 nanoseconds now used, 
the range biases for most GRARR passes would move to within one sigma of zero. 

The magnitude of the center frequency offset necessary to fully explain the 
GRARR range timing errors noted in this report lies well outside the operating re- 
gion of GEOS-A. However, it is felt that oscillator drift has had a significant 
affect on the range timing error figures. As an example of the errors caused by 
a frequency offset, consider a hypothetical shift S in the center frequency of the 
crystal oscillator aboard the spacecraft and assume no other change in shape of 
the transponder correction curve. 
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Let: 


S = 25 kHz « - 3300 meters/sec. 


The transponder delay correction A R T has the form: 

AR t = aR 2 +bR 

-8 -3 

where a = 7.18 x 10 , b = 0.332 x 10 and R is the measured range rate. 

Introducing a frequency shift S, the transponder delay correction becomes: 

AR = a(R + S) 2 + b(R + S) 
s 

Expanding and regrouping the above expression, AR becomes: 

s 

AR = AR„ + 2aSR + (aS 2 + bs) 

S X 

Hence the error in the transponder delay correction introduced by a frequency 
shift S will be: 

E = 2aSR + (aS 2 + bS) 
s 

Assuming S to be constant over the duration of a short arc pass, E takes 

s 

the form: 

E = AT R + B 
s 

where 

AT = 2 aS 
B - aS 2 + bS 

Using the values specified previously, the errors caused by neglecting the 
shift S would be: 

AT =-0.5 milliseconds 

B = - 0.3 meters 

Thus a timing error of - 0.5 milliseconds would occur if the transponder oscil- 
lator had shifted 25 kHz as shown by the dotted curve in Figure 1 and the total delay 
in the spacecraft transponder were still represented by the solid curve in Figure 1. 
As was stated previously, the range bias error caused by the shift is insignificant. 
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Another error model was postulated to explain range errors. It consisted of a 
term representing a servo lag error in the system and a timing error term and had 
the form: 

AR = K ? R + Kg R 

This model did not significantly reduce the RMS errors of the range residuals as 
shown in Figure 3 . 


RANGE RATE ERRORS 


The average range rate RMS value for GDAP for the three passes which re- 
ceived concentrated analysis was 2.6 centimeters per second. 

After the range rate residuals were corrected for ionospheric and tropospheric 
13 

refraction effects , the average RMS was 3.6 cm/sec. Since the average RMS in- 
creased after refraction correction, there was some uncertainty about the accuracy 
of the functional form used. It was known that there existed a 5% uncertainty in the 
tropospheric correction and a 25% uncertainty in the ionospheric correction. How- 
ever, after refraction correction, plots of the three sets of range rate residuals 
gave an "S" shaped curve similar to that of Figure 2. These results were essen- 
tially in agreement with those of other short arc orbital work using orbits defined 
by appropriately weighted SECOR, TRANET, GRARR and optical data 14 as well as 
with twenty-four hour long arc fits using GRARR data with optically defined orbits. 


Since there were some misgivings on the refraction correction, a regression 
analysis was done using the model: 


A R 


where 


C 


R 


= K 5 C R + K l® + K 2® 

_ E cos E 
2 

sin E 


where E is the elevation angle of the target and E is its elevation rate. 
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The R term, which represents a frequency error, was included since work done 

at D. Brown Associates (see reference 14) indicated that this term was significant. 

• • 

The R, representing a range rate timing error, was used since a consistent -0.2 
millisecond error was noted for all three passes using range rate error model (1) 
in Figure 3. 

The refraction parameter Kg which was to be determined was of the form: 


K 5 



where N was the total refractivity. 

The results of the regression analysis for the refraction error model coefficient, 
K 5 , showed a high correlation of 0. 89 between the range rate residuals and C R , in- 
dicating this term can fit the data and reduce the residuals. However, the values of 
K- were unrealistic in a physical sense and therefore it was concluded that it was 

D 

better to apply the theoretical refraction correction. 

The results for the R term, similarly indicated a high correlation of -.92 
between the range rate residuals and the range rate. Also, the use of a model of 
the form: 

AR = ^R 


reduced the average range rate residual RMS from 3.63 cm/sec to 1.36 cm/sec. 
In the physical GRARR system the coefficient can be derived as follows: 


Using 

(a) 


R 



where f D is twice the doppler frequency and f T is the up-link transmitter 
frequency. 

The representation of R shown above is an approximation sufficient to determine 
first order effects such as transmitter frequency offset. 


Differentiating (a) we get: 


AR 
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but substituting (a) into (b) 



hence: f 

T „ _ AC A T 
K 1 C " f T 

Since the same value of C, the speed of light, is used in both the Laser and the 

GRARR preprocessing, any error therein would cancel. Therefore would 

represent an error in the transmitter frequency. The carrier frequency error 

obtained from was approximately 1 part in 10 or 22.7 kc at a frequency of 

2271.9328 Me. The ground station oscillator which produces the carrier frequency 

7 

has an accuracy specification of 5 parts in 10 or approximately 1.1 kc for thirty 
days. However, drifts of as much as 4. 8 kc have been noted when the thirty day 
calibration has been performed. Therefore, it was concluded that if there were an 
error in the transmitter frequency or in its value used for preprocessing, it could 
not be of the magnitude suggested by and hence, other systematic effects must 
be present. Therefore, other error models for range rate were postulated. 

Since the curve in Figure 2 can be represented by many functional forms, ob- 
taining a satisfactory error model for range rate became a matter of trial and 
error. Reduction in the residual RMS, consistency of error model coefficients for 
the three passes and physical significance were the primary prerequisites for 
accepting an error model term. Since a -0.2 millisecond timing error appeared 
consistently in all previous error models using R, this term was used with a term 
representing a Type n servo lag. 

The model was of the form: 

(a) AR = K 3 *R*+ K 2 R 

where represents the timing error and Kg the servo lag error. 

The reduction in the average RMS for the three passes was from 3.63 cm/sec 
to 1.84 cm/sec and, again, a timing error of -0.2 milliseconds was observed. 
Since the reduction in RMS was not as good as was desired and the servo lag co- 
efficient was inconsistent, an R*R term was added to the model (a) in order to give 
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it the form of the differential of the transponder range delay curve, which had been 
tested previously, coupled with a servo lag term. 

AR = K 3 R’+ K 2 R + K 4 RR 

0 0 0 99 

where R is the servo lag model and K 2 R + K^R R is the differential form of the 
transponder range delay curve. The RR term was retained because of its signifi- 
cant contribution to the reduction of the R RMS. Its physical significance in terms 
of the range rate delay through the transponder has not yet been determined. 

The reduction of RMS from 3.63 cm/sec to 1.25 cm/sec was somewhat better 
than the previous model. However, it was noted theat the servo lag was the least 
significant term in the regression and reduced the RMS only .01 cm/sec. There- 
fore, a servo lag error was eliminated from further consideration in the range 
rate error modeling. 

The final range rate error model consisted of all terms which appeared signifi- 
cant in the previous models tested and was of the form: 

A R = K 5 C r + K^R + K 2 R + K 4 R R 

Using this model, a significant reduction in average RMS from 3.63 cm/sec to 
1.01 cm/sec was observed. 

As previously stated, using the Laser orbits as a standard, the accuracy of 
range rate could be determined to no better than 1 cm/sec. Therefore, this model 
seemed to define the range rate errors as well as could be expected. 

The most significant term in the regression was the RR term followed by the 
• * f 

R term, while the R and refraction (C D ) terms appeared to be of equal significance. 

n 

That is, when either the R term of the C R term was entered into the regression, the 
addition of the other term did not affect the fit at all. Therefore, both the model 

AR = k 5 c r + k 2 r + k 4 rr 

and the model 

AR = K 4 R + K 2 R + K 4 RR 

explained the range rate residual curve of Figure 2 to the expected accuracy. The 
coefficients K 2 and K 4 showed consistency, whereas the coefficients of C R and R did 
not. This makes it difficult to justify their presence in the error model, though the 
presence of one or the other contributed to a reduction in the RMS value . 
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ANTENNA MOTION DOPPLER 


A further correction to the range rate values due to motion of the antenna 
about the X-axis is suggested and will be implemented in further analysis done on 
the GRARR system. It was felt that the range rate accuracy attainable in this ex- 
periment did not warrant its inclusion; however, in work where better accuracy is 
required it will be added. The correction should be added to the range rate measure- 
ment and is given in meters per second by: 

AR = -1.17Y sin Y 

where Y is the Y angle of the GRARR antenna and Y is its rate of change. 

For GEOS-A the maximum value of the correction is about 0. 2 cm/ sec. 


ERROR IN PROPAGATION ANOMALY 

A tropospheric correction only is made to the range measurements in the GDAP. 
The functional form to be added to the range is given by: 


ar t = 


-2a 


2 

sin E + / sin E + 


4H„ 


where is the scale height of the troposphere, r Q is the earth radius, E is the 
elevation angle of the spacecraft and 

* ■ - <"o - l > [< n o - « r „ + H t] 
where n Q is the ground index of refraction. 

The nominal tropospheric zenith correction for range is approximately -2.0 + 0. 1 
meters. 


The range rate residuals are corrected for tropospheric and ionospheric 
refraction after being processed by GDAP. The index of refraction n Q used for the 
tropospheric refraction correction is calculated using the pressure, temperature 
and humidity at the station at the time of the pass. As was stated previously, this 
correction is known to be in error by as much as 5%. The ionospheric portion of 
the correction is known to be in error by as much as 25%. 
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The functional form for the total range rate refraction correction to be added to the 

13 

range rate is given by: 


AR 


E cos E 

. 2 „ 
sin E 





where N. is the ionospheric refractivity and E the elevation rate of the spacecraft. 
The terms containing and k give the tropospheric portion of the refraction cor- 
rection, k is a table look-up function of n . 15 

o 

A similar form 13 due to J, J. Freeman for the range ionospheric refraction 
correction is given by: 


AR = 


sin E 


/* N.dh 


cot 2 E 


/°°N.hdh 

J o i 


L o j 

The nominal ionospheric range refraction at zenith for passes at night is 0.5 + 
0.13 meters. 


There was some uncertainty about the refraction corrections used since the 
up-link carrier frequency of 2271. 9328 me from which doppler is derived is re- 
transmitted down-link with a frequency of 1705 me. The fact that the up-link 
carrier effectively traveling with phase velocity and the down-link subcarrier 
traveling with group velocity are at different frequencies had not been taken into 
account in the refraction correction. In order to remedy this situation in the future, 
equivalent frequencies for both range and range rate refraction corrections will be 
used rather than the up-link frequency now used. These frequencies are derived as 
follows. 


Both J^N.dh and J^N.hdh are functions of the maximum usable frequency at 

3000 km (MUF (3000)) , the plasma frequency, f Q , and the operating frequency of 
15 


The values of MUF (3000) and f Q for a given month at a given 


the system, f. 

position can be obtained from the CRPL Ionospheric Prediction Map 
the month and the position of the tracking system: 


Hence, given 
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/rv» = a (72-) and C N i Mh ■ b (t2 

where A and B are constant 

Therefore, given the radius of the earth, r , both the range and the range rate re- 
fraction corrections at a given elevation will be inversely proportional to the oper- 
ating frequency squared. 

Since the ranging sidetones travel with group velocity both up-link and down- 
link it is sufficient to average the up-link (AR^Jand down-link (AR d ) refraction 

corrections to get the equivalent refraction correction (AR ) for range. There- 

eq 

fore 


A R eq = — 5— (ARu + A Rd ) 

Using the up-link (f^)and down-link (f^) frequencies it is possible to calculate a 

single equivalent frequency (f ) to be used in the standard refractiQn correction 

eq 

model. 

Since for a given elevation and under given environmental conditions the re- 
fraction correction is inversely proportional to the operating frequency squared: 



For channel A, f^ = 2272 Me and f Q = 1705 Me 

Substituting in (a) 

f (range) = 1928 Me 
eq 

For GEOS-A at an elevation of 30 degrees the range refraction correction using the 
equivalent frequency would be approximately 0. 3 meters greater than that using the 
GRARR up-link frequency. 


- 24 - 



For range rate the derivation is similar however it is necessary to reverse the 

sign of the index of refraction used in the up-link refraction correction since the 

16 

carrier travels with phase velocity up-link. Thus: 

^ R EQ = ~ 2 ~ (^ R D ~ AR u) 

As in the range calculation: 



Substituting in (b) 

f (range rate) = 3648 Me 

For GEOS-A at an elevation of 30 degrees the range rate refraction correction 
using the equivalent frequency would be approximately 0. 1 cm/ sec less than that 
using the GRARR up-link frequency. 


RANGE RATE AVERAGING ERROR 


This correction was not used in the intercomparison investigation but will be 

17 

used in further processing of range rate data. This error is due to approximating 
the instantaneous value of range rate by the average value over an interval AT. 

For an orbit such as that of GEOS-A a maximum error of about 0. 7 cm/sec could 
occur in the range rate. The correction is added to the measured range rate and 
is given by: 


Ar 


R* ( A T) 2 
24 


where R is the third time derivative of range. 
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CONCLUSIONS AND SUMMARY 


The GRARR range bias of -5. 3 + 2. 5 meters and the range timing error of 
-2. 1 + 1. 2 milliseconds could be explained in part by inaccuracies in the transponder 
delay curve. More passes must be analyzed to determine the magnitude of this 
effect. 

A consistent range rate timing error of -0. 20 + .02 milliseconds was observed, 
but its cause has not been detected. The "S" shaped range rate residual curve 
observed in both short arc and long arc orbital work can be satisfactorily defined 
to 1 cm/ sec RMS by the model: 

AR = K,R + K„R + K.RR 
12 4 

This curve and, hence, the above model may be altered by further refinements in 
the range rate refraction correction. 

From this intercomparison experiment it appears that Laser orbits can be 
used to detect systematic errors in both the range and the range rate to about 2 
meters and 1 cm/sec respectively. After editing the Laser data, systematic effects 
therein seemed to be at a minimum. 

It is felt from the work done here that for future spacecrafts carrying trans- 
ponders, if a more thorough analysis is made of the delay and draft characteristics 
of the transponder, the GRARR system accuracy may exceed the 15 meter and 
10 Cm ^sec design goals for the range and range rate measurements. It would be 
helpful in determining transponder inaccuracies if the drift in the down-link fre- 
quency could be measured and recorded periodically. 
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TRANSPONDER DELAY (NANOSECONDS) 
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Figure 1. Range Error (Transponder Delay) vs. Doppler (Range Rate) 
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l RANGE ERROR MODELS Ij 




INITIAL RANGE RMS 6. 2 METERS 



MODEL 

FINAL RMS 

(1) 

A R 

= 

V 2 + K 6 R 

6.1 

(2) 

A R 

- 

K ? R + KgR 

6. 1 


RANGE RATE ERROR MODELS | 
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MODEL 

FINAL RMS 

(1) 

A R 

= 

2K 4 R R + K 2 R 

1.26 

(2) 

A R 

= 

K 5 C R + K 1 R + K 2 R 

1. 09 

(3) 

A R 

= 

K X R 

1. 36 

(4) 

A R 

- 

k 3 r’ + k 2 r 

1.84 

(5) 

A R 

= 

k 3 r*+ k 2 r + k 4 rr 

1. 25 

(6) 

A R 

= 

K & C r + K^R + K 2 R + K 4 R r 

1.01 

(7) 

A R 

= 

I^R + K 2 R + K 4 R R 

1.01 
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Figure 3. GRARR Error Models 
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Table 9 

Coefficients of GRARR Error Models 





= transponder delay term (in seconds /meter) 

= range rate refraction error (in centimeters /second) 
= range timing error (in milliseconds) 

= range servo lag error (in seconds) 

= transponder delay term (in seconds /meter) 
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